Binary projective measurement via linear optics and photon counting 
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We investigate the implementation of binary projective measurements with linear optics. This 
problem can be viewed as a single-shot discrimination of two orthogonal pure quantum states. We 
show that any two orthogonal states can be perfectly discriminated using only linear optics, photon 
counting, coherent ancillary states, and feedforward. The statement holds in the asymptotic limit 
of large number of these physical resources. 
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Projection measurements play an essential role in pho- 
tonic quantum-information protocols. In these applica- 
tions, generally, a projection onto superposition states or 
entangled states of optical fields is required. Physically, 
it is a highly nontrivial problem how to implement such 
a measurement. 

One plausible approach is to use linear optics and clas- 
sical feedforward associated with a partial measurement. 
For example, a universal quantum computation scheme 
for photonic-qubit states has been proposed, which uti- 
lizes only linear optics, photon counting, and highly en- 
tangled auxiliary states of n photons generated by prob- 
abilistic gate operations pj. In principle, it works with 
unit success probability in the asymptotic limit of large n. 
It is, however, still a nontrivial question how to prepare 
entangled ancillae even for modest n. 

In this paper, we discuss the linear optics implemen- 
tation of a measurement which effects a projection onto 
two orthogonal states {|\&), |$)}- This is equivalent to 
the problem of discriminating two orthogonal quantum 
signals {1*),]$)} unambiguously 0, II]. We show that, 
in the asymptotic limit of a large number of partial mea- 
surements, one can perfectly discriminate the two states 
with linear optics, photon counting, and feedforward, but 
without any non-classical auxiliary states. Even in the 
worst case, the average error probability of discrimination 
approaches zero with the scaling factor of TV -1 / 3 where 
N is the number of the partial measurements. Note that 
the signal space is two-dimensional but |\&) and |$) can 
be any physical states defined in a larger space, e.g. qubit 
states, continuous variable states, etc. 

Before discussing a linear optics implementation, it is 
worth mentioning a result concerning the distinguisha- 
bility of two orthogonal multi-partite states via local op- 
erations and classical communication (LOCC). The nec- 
essary condition for exact local distinguishability is that, 
after doing a measurement at some local site, every pos- 
sible remainingstates must be orthogonal to each other. 
Walgate et al. ,J] showed that there always exists a local 



projective measurement satisfying this orthogonality con- 
dition for any set of two orthogonal states. Thus one can 
perfectly discriminate them via a series of local projective 
measurements where the choice of the measurement basis 
at each local site is conditioned on the previous measure- 
ment outcomes. This result means that if one can show 
a physical scheme that can exactly discriminate any two 
orthogonal single-mode states, its sequential application 
can achieve an exact discrimination of any two orthogonal 
multi-mode states. In the following, therefore, we concen- 
trate on a discrimination of two single-mode states. 

An arbitrary set of two orthogonal single-mode states 
are described by 
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where \m) is an m-photon number state and (^ , |$) = 
J2m=o c mdm = 0. Figure ^ is the schematic of the mea- 
surement apparatus. The states are equally split into N 
modes by TV — 1 asymmetric beamsplitters Q , 

Bn-i,o(9n-i)Bn-2,o(8n-2) ■ ■ ■ Bi^ (8i)\0)® N 1 | v I')o 

= N BS \y) , 
where B,,o(^) = exp[6>;(aJao - a* 



(2) 

;aj)] 6] and tan0* = 
1/y/N — i. The input is symmetrically split to N modes 
with the effective power reflectance of 1 /N . Then, at each 
output port, one makes some measurement by using lin- 
ear optics and photon counters, where the information 
about the measurement outcome is fed forward to design 
the next measurement. It should be noted that this is a 
generalized version of the scheme so-called "Dolinar re- 
ceiver" 0, H El which was originally proposed as a physi- 
cal model attaining the minimum error discrimination of 
the binary coherent signals {|a), | — a)}. 

We briefly sketch how two states are discriminated 
by such a scheme in the limit of N — > oo and then 
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provide a rigorous proof. Suppose one inserts \^) or 
|$) into the first beamsplitter. For sufficiently small 
1/N, the reflectance of multi-photons can be neglected. 
The states after beamsplitting are approximated to 
be S 1 ,o(^i)|0) 1 |*)o « |0)i|%)o + iV- 1/2 |f)i|?7i)o, and 
Si,o(0i)|O)i|$) o « |0)i|^ )o + iV- 1/2 |f)i|^i)o, where, 
(770 1 1^0) + (vi\ v i)/N ~ 0, since a beamsplitting opera- 
tion is unitary. Then mode 1 is measured. The mea- 
surement here is required to maintain the orthogonality 
of any conditional outputs of |^) and |$). The local 
measurement satisfying this condition is described by a 
two-dimensional projective measurement, 
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ko> = AT p0 ||0> 

= A^ {|0)-(A + O(A 2 ))|l)}, (3) 
= M pl {(X*+O(X*))\0) + \l)}. (4) 
where, Apo and A/" p i are the normalization factors and 

^(.\{m\^)\ 2 -\{mM\ 3 ) 

Here, we have assumed |(^o|^i)| 2 — K^il^o)! 2 7^ which 
implies X oc l/\fN and thus we can take \X\ <C 1 in 
the limit of large N. The other case, i.e. K^ol^i)] 2 — 
I (?7x I ^0) 1 2 = 0, will be discussed later. Under this as- 
sumption, the projective measurement of Eqs. © and 
(0J can be implemented by the displacement operation 
D(/3i/y/N) and photon counting as shown in Fig. ^b). 
Since both the signal and displacement are sufficiently 
weak, the corresponding measurement vectors are de- 
scribed by 
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which can be same as Eqs. (0 and l@J by choosing ap- 
propriate (3\. 

The conditional states after the first measurement can 



be rewritten again as \^') = YTr, 



\m) and |<3>') 



J2m=o c ^m\ m )- Since Nbs splits a state symmetrically, 
one can repeat the same procedure for the remaining 
state with the second beamsplitter, the displacement op- 
eration D^fo/VN), where /?2 is conditioned on the previ- 
ous measurement outcome, and a photon counter. After 
repeating the same procedure to modes 1 to N — 1 with 
appropriate /Vs, the final states at mode contain with 
dominating weight at most one photon and are still or- 
thogonal to each other. As a consequence, applying the 
final (iV-th) displacement and photon counting, one can 
exactly discriminate |^) and |$) with unit success prob- 
ability. 

Now, we discuss the scheme rigorously, i.e. include the 
effects due to the multi-photon reflections at each beam- 
splitter, which contribute to the failure of the measure- 
ment or giving the incorrect decisions. Here, the input 
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FIG. 1: (a) JV-splitter, and (b) a measurement apparatus at 
each step. A displacement operation D(f3i / \/~N) is realized 
by combining the signal with a coherent state local oscillator 
\/3i/VN sin 9) via a beamsplitter with sufficiently small power 
reflectance of sin 2 0. 



states |\E') and |<£>) are always physical, that is, the aver- 
age power of them are finite. Moreover, we assume that 
the probability distribution in photon number of those 
states decreases exponentially as c m = c m e~ mx < 2 where 
x is a real positive number. The prior probabilities can 
be set to be equal without loss of generality. Finally we 
assume that the average powers of local oscillators always 
satisfy \(3i\ 2 < \C^\ 2 +0(1/N) where is a complex 
constant independent of N. 

After finishing a whole process of N measurement 
steps, one can classify the results according to the se- 
quential patterns of detected photon numbers. Let us 
denote the events in which all the photon counters de- 
tect zero or one photon by 'success' events and the others 
by 'failure' events. Because of the symmetry of the N- 
beamsplitting, the probability of detecting k photons at 
the i-th measurement on average over all possible mea 
surement patterns is given by 
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<a**/N h + 0(l/N k+1 ), 
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where l^,) = l)(C ( g i )| l I'), whose probability distri- 
bution still decreases exponentially in number basis 
(see Appendix A), and C™ ax is the maximum value of 

(* (3 j4 fc "ol*/3,)/ fc! for a11 * and possible inputs pU- The 
probability of resulting the failure event Pf a u is then 
bounded as 
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fail 



< (C£ iax 



/N 2 + 0(1/N 3 )) x N 
/N + 0{1/N 2 ), 



(9) 



which implies that Pf a u approaches to zero in the limit 
of large N, at least with the order of 1/N. 

Even if the detection is successful, the conditional 
states get slightly non-orthogonal after each measure- 
ment step. To see this, we revisit the first beamsplitter 
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Bi t o(9i). Let us describe the states after beamsplitting 
such that the orthogonal and non-orthogonal parts are 
separated as 

B llO (0i)|O)|*) = |0)|ry„)+iV- 1 / 2 |l)K)+^ 1 |2)| ?72 ) + ... 

= |0)ho)+iV- 1 /2|i ) | ??l)+i v- 3 / 2 |l)| ?7r ) 
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Bi,„(0i)|O>|$) = |0)k)+A^ 1/2 |l)h)+A^ 3/2 |l)k) 

oc 

+ J2N- k / 2 \k)W k ), (11) 



fe=2 

where the first two terms exactly satisfy the orthog- 
onality (t]o\vq) + (t]i\pi)/N = and the last terms 
represent the multi-photon reflection terms. Here, 

\m) = E"= cm(i - i/A0 m/ V>, N-i/M) = 

E"= 1 c TO (m/iV)V2(i_i/ J v)(™- 1 )/2| ro _i) )i v-V2| r?1 ) = 
E"=i c m (l - (1 - l/NT^lm - 1), and N~ 3 / 2 \ Vr ) = 
(l^i) — hi)) (l^'s are also obtained by replacing 
c rn with d TO ). The terms \r] r ), \v r ) and that for multi- 
photon reflections, which have been neglected in the pre- 
vious discussion, cause the residual non-orthogonality. 
Note that the leading terms of all vectors \rjk)'s and |ffc)'s 
are independent of N. Denote the i-th measurement op- 
eration as 



(12) 
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Then the conditional outputs after detecting zero and 
one photons at the first measurement are given by 



4%) = M + ^kg)], (13) 

e{ 1] \*) = M {/?i|%> + \m) + ^hS)} . ( 14 ) 

respectively, where Ao and M\ are the normalization fac- 
tors and the third terms \rf^. ')'s (i = 0, 1) come from \r) r ) 
and \r)kYs for fc > 2, and the terms in Eqs. © and (0) 
whose order is higher than 1 /iV 1 / 2 . The same outputs are 
obtained for |<I>) by replacing \rj n ) with \v n ). The first two 
terms in Eqs. I|13fl and (|14f> can be exactly orthogonal to 
those of |$) by choosing Px/y/N = (1 - Vl + X 2 )e iu /X, 
where X is obtained by substituting 1 770) , |^o) and 
appearing in Eqs. i|Ttj|) and (|TT|l . into Eq. ©. Since 
X oc l/VN as mentioned above, this choice of Pi always 
satisfy the constraint on the average power of the local 
oscillator, |/3i| 2 < \Cp 1 1 2 + 0{1/N). However, we have to 
care of the fact that, in both events, the total conditional 
states in Eqs. and (|14|) are no longer orthogonal due 
to their third terms. 

Now, suppose that the same strategy is applied to the 
choice of P2 for the second measurement step. After the 
second measurement, the states are mapped into the new 
one with orthogonal and non-orthogonal terms, where the 
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FIG. 2: The original scheme (a) can be transformed into (b) 
where the total input photon number is the sum of those of 
two input states. 



latter has two parts, i.e. contributions from the first and 
second measurements. Note that the leading order of 
prefactors of \ry^ ) with respect to 1/N does not change 
during the measurement process, as also the leading fac- 
tors of does not change in the mapping in Eqs. H13|) 
and (JT2J. Eventually, after repeating N— 1 measurement 
steps in a similar way, if all the photon counters detected 
zero or one photons, one obtains the conditional output 
consists of the orthogonal term and N—l non-orthogonal 
terms stemmed from each measurement as 
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where the first term is exactly orthogonal to that of 
l*^ 1 )), while is the residual non-orthogonal term 

coming from \rj^ ). J^ -1 ) and J^" 1 ) are the numbers 
of the events of detecting zero and one photon, respec- 
tively, and thus J^- 1 ) + J( JV - 1 ) =N-1. 

Let us denote the final iV-th measurement by \Du) = 
(Pn/ '\N)\k) (k = 0, 1). Suppose that Pn is designed 
such that \Dq) and |£>i) are the same as the orthogonal 
terms in l^^ -1 )) and |<I > ( Ar ~ 1 )), respectively, up to the 
order of 1/N 1 / 2 (the higher order terms contribute to 
the detection error). Then the error probability P£Q = 
I (Dil*^- 1 )) | 2 is given by 
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where + jW = N. The leading order of (Di\H^) 
is independent of N for every j and k. 

One can estimate the order of JW by count- 
ing the total amount of photons put into the sys- 
tem since the number of the total photon is equal 
to that of detectors. Although photons are supplied 
by the input state and N displacements in the orig- 
inal configuration, one can simplify it into the one 
with only two inputs, D(/3o)|^) and the coherent state 
|/3aux), by adding some linear optics as illustrated in 
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Fig. Here, with the relation D A (a)D B (P)B AB (9) = 
Bab(&)Da{o: cos 9 — (3 sin 9) D B (a sin 9 + /3cos#), one 

finds \(3 \ 2 = | Eti A/^l 2 and |/3 aU x| 2 = Eli \Pi\ 2 /N- 
|/3 | 2 , where these are bounded as \(3 \ 2 = Co + 0(1/N) 
and |/3 a ux| 2 = C aux + 0(1/N) due to the constraint on 
|/?i| 2 's. C and C aux are constants independent of iV. 

The probability of having n photons in total is given by 
p (n) = El= P si g (n-m)P^(m) = C P e~ nx + 0(1/N). 
Here the photon number statistics of two inputs, P S ig( m ) 
and Paux(TO) are exponential and Poissonian, which eas- 
ily implies that P(n) decreases exponentially with re- 
sepect to n (see Appendix C). Therefore, one can bound 
JW by some constant Cj with exponentially small ex- 
ception as 



Prob J (A,) <Cj + 0(1/ N) + Ni 

> 1 - C P exp[-(Cj + 0(l/N) + 
= 1 - C P e- Nt + 0(1/N) 
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(17) 



where e can be arbitrarily small for large N . Eventually, 
substituting it and P N > < N into Eq. (|16|l . one obtains 



(Di\H ) 



(Dim 



N 2 ' N 

< C E /N 2 + 0(1/N 3 ) + eO(l/N) + e 2 , (18) 

where {D x \H k ) m = ^{D^H^/L (L = /W and jW 
for k = 0, 1, respectively), and Ce is some constant in- 
dependent of N . In a similar manner, the same bound 
is derived for = \{D \& N ~V)\ 2 . Then, summing 

over all detection patterns, the average error probability 
is bounded as 
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< C/N + 0(l/N 2 )+0(l/N)e + e 2 
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(19) 



where C is some constant and P(H) is the probability 
to observe the measurement sequence pattern jj. As a 
consequence, in the limit of N — ► oo, one can discriminate 
\^>) and |$) with unit probability. 

Finally, we discuss the case |(?7o|^i)| 2 — K^iko)! 2 = 
in Eq. (JSJ, in which the desirable local measurement 
can not be implemented by a displacement and pho- 
ton counting. Here, let us consider the projection mea- 
surement consisting of slightly perturbed vectors |IIo) = 
VT~5\y) - VS\^>) and |ni) = VT~6\<I>) + VS\V) with 
a perturbation parameter S. One can design such a 
measurement by the previous strategy with the total er- 
ror probability of P*°* = C/N 1 -^ + 0{1/N 2 - 3A ) + 
0(l/JV 1 - SA / a )e + 0(N A )e 2 , where A = - \og N 5. This 
device can discriminate the original states \^>) and |$) 



with the average error probability of 

p- = l-ii-p^mom' + mmi 2 )^ 

= d/N A + 0(1/N 2A ) + C 2 /N X - 2A + 0{l/N 2 - 3A ) 
+0(\/N 1 - 3A ' 2 )e + 0{N A )e 2 (20) 



In the asymptotic limit of large N , this is minimized 
with A = 1/3 and then we obtain P^ r = C/N 1 / 3 + 
0{1/N 2 ' 3 ) + 0{1/N l ' 2 )e + 0{N x / 3 )e 2 which still con- 
verges to zero. 

In summary, we have proved that arbitrary two or- 
thogonal pure states can be perfectly discriminated by 
linear optics tools without using any non-classical ancil- 
lary states in the asymptotic limit of N — ► oo where N 
is the number of the detections and feedforwards. It im- 
plies that, in principle, one can implement arbitrary pro- 
jection measurement in any two-dimensional signal space 
by these tools. The resources discussed here are mostly 
available with current technology. We also showed a con- 
crete designing strategy of a linear optics circuit to at- 
tain this bound for a given iV and thus it can be directly 
applied for various quantum information protocols that 
require binary projection measurements. The remaining 
question is whether one can apply a same approach to 
the problem of more than three states discrimination. 
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APPENDIX A: PHOTON NUMBER STATISTICS 
OF THE DISPLACED STATE 

In this appendix, we show that if the photon number 
distribution of the initial state is exponential, then that 
of its displaced state is also bounded by exponentially de- 
creasing function. For this purpose we use the following 
three formulae; 

(1) The number basis components of the dis- 
placement operator [11]; 



(n\Dtt)\m) = \ '—tr- m e-l«l 2 / 2 z4r™)(|£| 2 ),(Al) 
V n! 



for (n > m) and 
(n\D(0\m) = 



'^ ( _ r) m-n e -|eiVa £ (m-n) ( | f|2)j 

(A2) 



for (n < m), where L„\x) is the associated Laguerre 
polynomial defined by 
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where Lffl (x) = L n (x) is the Laguerre polynomial and 

_. 

dx l 



^jL n (x) = (-l)^(z). 



(A4) 



Proof. We basically follow the proof given in 0. 
To calculate {n\D(£)\m), it is helpful to see (n\D(£)\n), 
which is given by 

(n\D(0\n) = (n\ exp (£a f - \n) 
= e-l«l 2 / 2 (n|e^ at e-«* a |n) 
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Then we obtain 
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Therefore, replacing n — I with m in Eq. (HI with 
Eq. <|A4|> . we can derive Eq. (|A1|) . 

(2) Bound on the associated Laguerre polyno- 
mials |l2[: 
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where a; > and a is an integer. 

Proof. From Eq. (|A5|I . the absolute value of the La- 
guerre polynomial L n {x) with x > is bounded by 



\L n {x)\ = e x ' 2 (nlDix 1 / 2 )^) 
< e x/2 . 



(A8) 



To extend it to the associated Laguerre polynomial, we 
use the relation 

L^{x)=j^( a + k - l ')L^ k (x) 1 (A9) 
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which can be derived as 



n—k 



a + k — l\ ^— ^ (n — k\ (—x) 1 



1=0 



k=0 
n n—l 

EE 

2=0 k=0 
n / \ I \l 

y-v ( n + a \ {—x) 



I J l\ 



a + k-l\/n-k\ (-x) 1 
a-1 )( I ) 11 



1=0 

where the formula 

n—l 



n-lj l\ 



E 

k=0 



a + k\ (n—k 



I 



n + a + 1 
n — I 



(A10) 



(All) 



has been utilized. Eventually, Eqs. ljA8j) . (|X9|) and (|A11|) 

imply 
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(3) Inequality for the binomial distribution [13] ; 

I y v (l - y) n ~ v < exp \-2n(y - vjnf\ , (A13) 



where n > v and < y < 1. 
Proof. Define q = v jn and 



/(!,) = y u (l - y )«- e -2»(- 9 ) 2 
F( y ) - n-Hnfix). 



(A14) 



Then 



= qlny+(l-q)ln{l-x) + 2(x-q) 2 , 
{q-y){l-2yf 



F'{y) 



2/(1 - 2/) 



(A15) 



and thus F(y) takes its maximum at y = q. Also, the 
same for f(y). Therefore, f(y) < /(g), i.e. 

y v (l - y) n -" e 2n ^-^ 2 < q v (l - q) n - v , (A16) 
and thus 

("W - y) n - v e^-tf < - g)"-" < 1, 

(A17) 

which completes the proof. 
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Derivation of the displaced state. Now we derive the main statement of this appendix. We assume that 
can be written as 



!*> = £ 



c m e- mx ' 2 \m) 



(A18) 
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Now, let us calculate (n\D(f3)\^). 
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1 \n J (m — n)! 



(A19) 



The last line follows from the Cauchy-Schwarz inequality. 

Introducing a real parameter q which satisfies e~ x < q < 1, the first term of Eq. (|A19|) is then bounded as 



^e- nx ' 2 \c max \ 



E 



,m=0 



n\(l/3| 2 e x ) 



2 „ a: \7i— m 



1/2 



\to J (n — to)! 



v^e— / 2 |c max | E '/:" 



g|/?| 2 e^" "' 



< 

q 



-x \ ™/2 



(n — to)! \ 1 — q 

^rw-.r-exp 



1/2 



K 171 — 



g|/3|V 
1-9 



1/2 



-x\ ™/2 



|c max | exp 



q\(3\ 2 e x 



2(1 -q) 



(A20) 



and thus it decreases exponentially as n increases. Also, for the second term, one obtains 

1/2 



Ke- nx ' 2 \c^\ 



E 



K m— n+ 
( oo 



2 „ — x\m—n 



l e~ x ) 



^e- nx / 2 \c ma J 



1 \n J (to — n)\ 



^ m— n+1 



(to — n)! \ 1 — <? 



1/2 



< x/n 



-x \ ™/2 



v m— n+1 



E n<r(i-9r-*«p 



«|/3|V 



1/2 



^/2 



|Cmax| exp 



9l/?| 2 e- 



2(1-9) 



1-9 

E ex p _2m (' 7 "^) 



v m— n+1 



1/2 



(A21) 
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Since the last exponential term decreases exponentially as m increase at least in the limit of m S> n, the sum always 
converges within a finite value, which means that Eq. (|A21|I itself also decreases exponentially as n increases. As a 
consequence, these results imply that Eq. (|A19(I decreases exponentially as n increases. 



APPENDIX B: DERIVATION OF THE INEQUALITY © 

In this appendix, we derive the inequality JSJ). 



A(A/^)£jv-i,o(0iv-i) • • • A,o(0i) • • • Si :O (^i)|0)^ A '- 1 |*)o | 
A(A/^)Siv-i,o(^iv-i) • • • 4+1,0 (0i+i)4-i,o(0i) • • • B 1 ,o(0 2 )S l ,o(^i)|O)^ Ar ^ 1 |*)o 
A(/3 l /ViV)B 4 ,o(^i)|0) l |*)o| 2 

e -|/3i| 2 /2Af g/ 3 i a i t /\/lV e -/3*a i /v / lV e -a i t ao taring- In cos 0i (at a< -a+a ) gSiaJ tan0i IqNJ^Nq 



e -|ft| 2 /2A r e /3ial/v / ]V e -a i t a /v / lV^T e /3* a /y/N(N-l) e -/3? ai/v^Vgln^/l-l/JVa+ao . |$^ ( 
1 



-|/3i| 2 /2JV 



E 

J=0 



fc! 



ft 



Of) 



-|ft| 2 /2JV 



2 

[3i - a 



e 0;a o /y/N(N-l) e ln^l-l/Nala o iQ\ 



e a+a ln^/l-l/JV ( M» ^_»0 ] g/3*a /JV|^ 



(Bl) 



where cos6»i = ^1 - 1/N, sin6»i = 1/-//V, |ft| 2 < 



the relation e aS 'e' 3 ^ 



' ft | 2 + 0(l/iV) and = D(C A )|*). We have used 

e QOj e a/3 from line 5 to 6, and 

jc- - - - from line 7 to 8, where a, (3, 

and 4> are complex numbers. These relations are directly 
obtained from the commutation relation [aj,aj] = 1. 

The remaining task is to show that (vf^ |aJ k a k \^^ i )/k\ 
is always finite, i.e. 



tfc-fc 



ft; 



N fc fc! 



< 



(B2) 



with a constant (7™ ax 



. Here, we replace j3i by fi for sim- 
plicity. As shown in Appendix A, the photon number 
distribution of j^) decreases exponentially. Denote 



OO 

E' 

m=0 



-mx/2 I 



(B3) 



The absolute of complex coefficients 6 m 's are always in 
between and some constant due to the normalization 
constraint and let us denote the constant as |6 max |. Then 



has 



a fe |*/3> = E ~ hm€ 



c/2 



m—k 

oo 



(m — fc)! 



1 771 — fc) 



m=0 



and thus 



Aftfc! 2^ 1° 



N fc fc! 



2 (m_+fe)! a 

m+fe I . 



m! 



< 



I?) | 2 p - fe;t 
I "max | c 



fc! 



A^ fe fc! (l-e" x ) 



1 |6„ 



fc+i 



iVfc i 



1 - e~ x 



(B5) 



This bound depends on 6 max and x i.e. the state l^ft). 
Therefore, maximizing the rhs of this inequality for all 
\^Bi) and denoting the maximum value as C™ ax /A^ fe , one 
obtains Eq. (|B2|) . 
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APPENDIX C: TOTAL PHOTON NUMBER which decreases exponentially as n increases. 

STATISTICS 

The exponential and Poissonian distributions are de- 
scribed as 

P E {m) = C E e- mx , (CI) 

and 

P P (m) = ^re~ Cp , (C2) 
ml 

respectively. Then the distribution of the total photon 
number is given by 

n 



m— 



e C P y ^ CEe -(n-m) X 

' ml 

m=0 



(C P e* 



ml 

m=0 



< C E e Cp ^-^e- nx , (C3) 
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